Pseudospectral Solution of Near-Singular Problems using Numerical Coordinate Transformations Based on Adaptivity

The work presented here describes a method of coordinate transformation that enables spectral methods to be applied efficiently to differential problems with steep solutions. The approach makes use of the adaptive finite difference method presented by Huang and Sloan [SIAM J. Sci. Comput., 15 (1994), pp. 776–797]. This method is applied on a coarse grid to obtain a rough approximation of the solution and a suitable adapted mesh. The adaptive finite difference solution permits the construction of a smooth coordinate transformation that relates the computational space to the physical space. The map between the spaces is based on Chebyshev polynomial interpolation. Finally, the standard pseudospectral (PS) method is applied to the transformed differential problem to obtain highly accurate, nonoscillatory numerical solutions. Numerical results are presented for steady problems in one and two space dimensions.